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ABSTRACT 

We present a novel method for interpreting observations of high-mass X-ray binaries 
(HMXBs) based on a combination of spectroscopic data and numerical results from a ra¬ 
diation hydrodynamic model of stellar winds. We calculate synthetic Doppler tomograms of 
predicted emission in low/hard and high/soft X-ray states and compare them with Doppler 
tomograms produced using spectra of Cygnus X-1, a prototype of HMXBs. Emission from 
HMXBs is determined by local conditions within the circumstellar medium, namely density, 
temperature, and ionization state. These quantities depend strongly on the X-ray state of the 
systems. By increasing intensity of an X-ray emission produced by the compact companion 
in the HMXB model, we achieved a complete redistribution of the circumstellar medium in 
the vicinity of the modelled system. These changes (which simulate the transitions between 
two major spectral states) are also apparent in the synthetic Doppler tomograms which are in 
good agreement with the observations. 

Key words: methods: numerical - techniques: spectroscopic - stars: mass-loss - circumstel¬ 
lar matter. 


1 INTRODUCTION 

Cygnus X-1 (VI357 Cyg; HDE 226868) is one of the brightest and 
widely observed X-ray sources in the sky, and an archetype of high- 
mass X-ray binaries (HMXBs). The optical component - an 09.7 
lab supergiant (Gies & Bolton 1986) - and the compact component 
- most likely to be a stellar mass black hole - revolve around the 
centre of mass of the system with orbital period 5.6 d. The intense 
X-ray emission from the compact companion is thought to be pow¬ 
ered by the accretion of the circumstellar matter originating from 
the supergiant. A further period of 294 d (Priedhorsky et al. 1983), 
discovered in both optical and X-ray data, is thought to be associ¬ 
ated with precession of the accretion disc. The system undergoes 
aperiodic transitions among several major X-ray states represented 
by a global flux and spectral change. The first transition was ob¬ 
served by Tananbaum et al. (1972): the soft X-ray flux (2-6 keV) 
decreased by a factor of 4, while the hard flux (10-20 keV) in¬ 
creased by a factor of 2, and radio emission turned on. Since the 
X-ray emission is produced through the accretion process of the 
stellar wind of the donor star on to the compact component it is 
safe to assume that such a transition is associated with a global re¬ 
distribution of the circumstellar matter in the system. 

Despite precipitous advances in observational techniques and 
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instruments, the angular resolution of even the largest telescopes 
is still insufficient to observe the accretion process in HMXBs 
through direct imaging. We thus have to employ indirect methods 
such as Doppler tomography which can resolve accretion flows in 
binaries on micro-arcsecond scales by using time-resolved spec¬ 
troscopy. This technique utilizes binary spectral line profiles taken 
at a series of different orbital phases and transforms them into a 
distribution of emission over the binary. 

In Sec. 2, we give a detailed description of optical spectra that 
we have available, and the method we used to determine the X-ray 
state of the system for individual observations. Spectra were ob¬ 
tained in each of the major X-ray states of Cygnus X-1. We analyse 
the states separately in order to determine any structural changes of 
the stellar wind during transitions. In Sec. 3, we compare results ob¬ 
tained from the Doppler analysis of the spectra compartmentalized 
according to the X-ray state. We also discuss influence of length of 
the spectra-gathering period on the quality of the resulting tomo¬ 
grams. Sec. 4 describes input parameters of the radiation hydrody¬ 
namic model (Cechura & Hadrava 2015, referred to below as CH) 
of the stellar wind in Cygnus X-1. We use the predicted distribution 
of the circumstellar matter to create synthetic tomograms for indi¬ 
vidual X-ray states, and compare them with tomograms produced 
using H-alpha line. 


2 OBSERVATIONS 

In order to perform the Doppler mapping of Cygnus X-1, we have 
collected spectroscopic data from the Perek 2.05-m telescope at 
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Ondrejov observatory of the Astronomical Institute of the Academy 
of Sciences of the Czech Republic. We also included the spec¬ 
tra obtained at 2.16-m telescope at Xinglong station of National 
Astronomical Observatories, China (NAOC; Yan et al. 2008). The 
data consist of isolated measurements as well as several series of 
observations, and cover a period between 2003 and 2013. At the 
Ondrejov observatory, 89 spectra were obtained with CCD in the 
focal point of 700 mm camera of Coude spectrograph. A typical 
resolution is about 0.25 A per pixel. At the Xinglong station, an ad¬ 
ditional 18 spectra were acquired with a CCD grating spectrograph 
at Cassegrain focal point of the telescope with an intermediate res¬ 
olution of 1.22 A per pixel. The spectra were obtained during the 
low/hard, high/soft and intermediate state of Cygnus X-1. The com¬ 
plete list of the spectra used in the Doppler analysis can be found 
in appendix, including observational date, UT time of the begin¬ 
ning of the observation, exposure time, heliocentric Julian date of 
the centre of the observation interval, wavelength range, and spec¬ 
tral state. Orbital phase is also indicated, the ephemeris of the infe¬ 
rior coniunction of the companion star is adopted from Gies et al. 
(2003), 

2451730.449 + 5.599829^ . 

The observations were reduced using standard reduction rou¬ 
tines in IRAF. All spectra were bias-subtracted, flat-field corrected, 
and had cosmic rays removed. After cleaning and wavelength cal¬ 
ibration, the spectra were converted to ASCII format to be further 
processed in MOLLY, a one-dimensional spectrum analysis pro¬ 
gram developed by Tom Marsh. In MOLLY, the spectra were cor¬ 
rected for the heliocentric effects, and normalized by fitting poly¬ 
nomials by which the corresponding spectra were divided. 

The most prominent feature in the observed spectral range is a 
variable Ha line profile that is strongly dependent on orbital phase 
and spectral state. It is apparent from the observations that the spec¬ 
tra taken in the low/hard state have a strong emission in the whole 
Ha line profile, while in the high/soft state period, the emission 
remains only in the long-wavelength wing of the line and the short- 
wavelength side of the line profiles reveals an absorption, as is typ¬ 
ical for the P-Cygni profiles of stars with intense mass-loss via a 
stellar wind. 

The He I A6678 line is practically free of emission in both 
states. It means that this line may enable us to reliably measure 
radial velocities (RVs) of the primary component to get a constraint 
on the orbital parameters of the system. 

One of the pronounced spectral features is the absorption at 
A = 6614 A. It is actually an example of a diffuse interstellar band 
(DIB) in our spectra. The carrier molecule of these interstellar ab¬ 
sorption bands is still unknown (Herbig 1995; Krelowski 2003). 
The identification of the carriers of these bands remains an impor¬ 
tant problem in astronomy to date, and the current consensus on 
the nature of the carriers is that they are probably carbon-bearing 
molecules that reside ubiquitously in the interstellar gas (Ehren- 
freund & Charnley 2000). The most promising carrier candidates 
are carbon chains, polycyclic aromatic hydrocarbons (PAHs), and 
fullerenes (Salama et al. 1996, 1999; Foing & Ehrenfreund 1994; 
Schulz et al. 2000; Motylewski et al. 2000). The DIB at A = 6614 
A has a finer structure consisting of three absorption lines (Galazut- 
dinov et al. 2002) which are, however, below the resolving power 
of our instruments. Since DIBs do not follow the orbital motion of 
the binary, we use them to increase accuracy in wavelength scaling. 
For this purpose, we calculated RVs of the systemic shifts for each 
spectra using KOREL, a code employing a method of Fourier dis¬ 


entangling of spectra of binary and multiple stars (cf. e.g. Hadrava 
2009, 2012; Yan et al. 2008). 

The spectra are of good quality in terms of spectral resolu¬ 
tion which makes them suitable for Doppler tomography. However, 
in the terms of time coverage, things get more complicated. It is 
assumed in Doppler tomography that the distribution of emission 
from the system is constant in the course of the whole observational 
campaign. However, Cygnus X-1 varies also on time-scale shorter 
than its orbital period Porb — 5.6 d. The characteristic dynamical 
time-scale of the circumstellar matter is of the order of Porb, which 
is also the minimum duration of observational run needed to cover 
all orbital phases. Because of this difficulty we thus suppose that, in 
spite of the short-time variations caused most likely by a turbulent 
motion of gas, the fundamental distribution of emission remains 
stable and unique for a particular spectral state. This assumption 
allows us to divide all observations in groups according to their 
spectral state. With two exceptions of observational series in 2003 
and 2013, most of our data are scattered throughout a period of 
almost a decade. These observations spread over many temporary 
fluctuations provide us with a sufficient phase coverage to perform 
the Doppler tomography and to find the mean distributions in the 
low/hard and high/soft states. In the case of the observational series 
in 2003 and 2013, we gathered, in the period of several months, 
enough spectra for the independent Doppler analysis. Even though 
the orbital phase coverage is not optimal, we can still make some 
general conclusions and compare the results with the tomograms 
obtained from the whole period 2003-2013. 

2.1 Determination of X-ray state of Cygnus X-1 

To determine the X-ray spectral state of Cygnus X-1 we need to 
use X-ray observational data. In the past decade, state information 
was readily available using the All Sky Monitor on the Rossi X- 
ray Timing Explorer (Levine et al. 1996, RXTE-ASM) and regu¬ 
lar pointed monitoring observations. Unfortunately, RXTE ceased 
science operation on 2012 January 5 (MJD 55931) but even prior 
to this date since January 2010, ASM experienced an instrumental 
decline which made the observed data unsuitable for determina¬ 
tion of the spectral state. We therefore had to combine X-ray data 
from multiple sources. Here, various state definitions exist, which 
use measured count rates or colour (Remillard 2005; Gies et al. 
2008), or sophisticated mapping between these measurements and 
spectral parameters (cf. e.g. Ibragimov et al. 2007; Zdziarski et al. 
2011). The former prescription is easy to use, but is very instru¬ 
ment specific and cannot be easily translated to other X-ray all sky 
monitors. The latter approach requires a sophisticated knowledge 
of the instrumentation of all sky monitors, as well as of the detailed 
spectral modelling. Furthermore, the previously used state defini¬ 
tions are slightly inconsistent among themselves. The source be¬ 
haviour from early 1996 (RXTE launch) until the end of 2004 has 
been discussed by Wilms et al. (2006), who used a crude definition 
of ASM-based states using only the count rate. This classification 
is sufficient to distinguish main activity patterns, and also consis¬ 
tent with more detailed studies (Zdziarski et al. 2011). But in this 
work, we use a novel approach to classify states of Cygnus X-1 us¬ 
ing all sky monitors RXTE-ASM, MAXI, and Swift-BAT based on 
16 years of pointed RXTE observations and suggested by Grinberg 
et al. (2013). 

MAXI is an all-sky monitor on board the Japanese module of 
the International Space Station (Matsuoka et al. 2009). Light curves 
from the Gas Slit Camera detector (GSC) are available in three en¬ 
ergy bands (2-4 keV, 4-10 keV, and 10-20 keV) on a dedicated 
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State 

ASM-based 


MAXI-based BAT-based 

Low/hard 

20Vc^ 55{h-ho) 


CM < 1.4/lM 

Intermediate 

c > 20 A 55(/i — ho) < c ^ 350(/i - 

- ho) 

1.4/iM < Cm ^ 8/3/im 

High/soft 

c > 20 A c > 350(/i — ho) 


8/3/im < Cm Cb ^ 0.09 


Table 1. X-ray state definitions according to different all sky monitors. ASM 1.5-12 keV count rate c is in counts , and h is ASM hardness defined as (5-12 
keV/1.5-3 keV), ho = 0.28. MAXI 2-4 keV count rate cm is in counts and Hm is MAXI hardness defined as (4-10 keV/2-4 keV). BAT normalized 
15-50 keV count rate cb is in counts cm“^ s“^. Distinguishing between hard and intermediate state is not possible from the BAT light curves alone. The 
source is defined as in BAT-based hard or intermediate state for cb > 0.09 counts cm“^ (Grinberg et al. 2013). 


website^. MAXI light curves show prolonged gaps of several days 
due to observational constraints. MAXI coverage started on MJD 
55058 (2009 August 15). 

Swift-BAT is sensitive in the 15-150 keV regime (Barthelmy 
et al. 2005; Krimm et al. 2013). Satellite-orbit averaged light curves 
in the 15-50 keV energy band from this coded mask instrument are 
available on a dedicated website^. BAT coverage started on MJD 
53414 (2005 mid-February). 

Table 1 shows the determination of state based on the data 
from different all-sky monitors taken from (Grinberg et al. 2013). 
A graphical representation associating light-curve behaviour with 
spectral state of Cygnus X-1 is given in Fig. 1. The hrst two panels 
show the RXTE-ASM data coloured aceording to their established 
spectral states (blue is for the low/hard, green for the intermediate 
and red for the high/soft state). The hrst panel represents the sum of 
counts in ASM 1.5-12 keV energy range; the seeond panel shows 
the hardness of spectra. The hardness is dehned as a ratio of count 
rates in C band (5.0-12 keV) and A band (1.5-3.0 keV). Due to 
the deterioration of the ASM following MJD 55200 (depicted by 
the dotted line) no data from the ASM after this date (grey colour) 
is used for our analysis. The third and fourth panels show the data 
from MAXI energy band 2-20 keV and Swift-BAT energy band 
15-50 keV, respectively. Black represents the low/hard and inter¬ 
mediate state which are indistinguishable from the Swift-BAT data 
alone. The lowest panel displays the distribution of optical spectra 
available for our analysis. The Ondfejov spectra he below while the 
Xinglong data from 2004 and 2006 are shifted upward. The data are 
scattered over a period since 2003 to 2013. The best series we have 
available are 2003 with 35 spectra and 2013 with 20 spectra. The 
individual observations are eoloured aceording to their correspond¬ 
ing spectral states. 

In Fig. 2, we show distribution of optical spectra from 
Ondfejov and Xinglong in orbital phase of Cygnus X-1 with re¬ 
spect to their eorresponding X-ray states. We assigned 44 spectra 
to the low/hard state (upper panel), 19 spectra to the intermediate 
state (middle panel), and 44 remaining spectra to the high/soft state 
(lower panel). Unfortunately, due to the cadence of our observa¬ 
tions, the effective phase coverage is worse than a truly random 
sample. Despite the clustering problem, the orbital coverage is suf- 
hcient to use for Doppler analysis in the low/hard and high/soft 
states. We highlight (longer lines) 17 spectra taken in 2003 during 
the low-hard state and 17 spectra taken in 2013 corresponding to 
the high-soft state. 

Fig. 3 shows the evolution of the main spectral features of 
HDE 2266868 over an orbital period. The left panel shows se¬ 
lected optical spectra from 2003 that correspond to the low/hard 
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state of Cygnus X-1. The corresponding observation dates and or¬ 
bital phases calculated according to the above given ephemeris are 
written on the left and right side of each spectrum, respectively. 
Upon close examination, the trailed spectra in Ka were found to 
consist of two superimposed components, shifted in orbital phase. 
The hrst component corresponds to a P-Cygni prohle of an expand¬ 
ing stellar wind coming from the donor star. The second component 
is consistent with a presence of strongly variable circumstellar mat¬ 
ter in the system. The spectra obtained in 2013, all consistent with 
the high/soft state, are shown on the right panel. In the high/soft 
state, only P-Cygni component varying with the orbital movement 
of the donor is detected. There are no traces of the companion or 
the circumstellar matter in the He I A6678 in either spectral state. 


3 DOPPLER MAPPING 

To analyse our optical spectra, we use MODMAP - a program for 
modulated Doppler tomography using MEM introduced and devel¬ 
oped by Steeghs (2003). A grid of pixels spanning velocity space 
is adjusted to achieve a target goodness of ht measured with 
while simultaneously maximizing the entropy of the image. A re- 
hned form of entropy which measures departure from a default im¬ 
age is used. Eor an excellent summary of the method, see Marsh 
( 2001 ). 

Eor Cygnus X-1, we used the systemic RV 7 = —2.7 km s“^, 
published by Gontcharov (2006). The systemic velocity needs to be 
known in order to calculate the correct RV curve of a line source 
with velocity (vx,Vy). If one attempts to reconstruct a data set with 
an incorrect 7 , it will be difficult to achieve a good ht to the data 
since MODMAP is not tracing the correct velocity curves (ef. Marsh 
& Horne 1988). 

We adopt the orbital period Porb = 5.599829 zb 0.000016 d 
calculated by Brocksopp et al. (1999). A high precision in the de¬ 
termination of orbital period is needed since we attempt to analyse 
the phase resolved data spanning over a decade. If the period is not 
sufficiently well determined, the tomography will fail completely. 
Less severe consequences would have an incorreet determination 
of zero-point ephemeris. In this case, the resulted tomogram would 
be inaccurately positioned. The mass donor star and the compact 
companion would be rotated by an unknown angle with respect to 
the emission. 


3.1 The low/hard state 

Here, we present the results of the Doppler tomography of Ka 
A6562.76 for the low/hard state of Cygnus X-1. Eirst, we analyse 
all the spectra (44) corresponding to the low/hard state. Then, we 
run the analysis on 2003 data only (17 spectra in total) to assess the 
influence of length of the spectra-gathering period on the quality 
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Figure 1. X-ray light curves (RXTE-ASM, MAXI, and Swift-BAI) of Cygnus X-1 in the highest available resolution. ASM hardness is calculated by dividing 
count rates in C band (5.0-12 keV) by count rates in A band (1.5-3.0 keV). Vertical dotted line represents starting point of instrumental decline of ASM. Blue, 
green and red colours represent X-ray spectral states of individual measurements classihed by respective classihcation for the given instruments as described 
in Table 1: blue represents the low/hard, green the intermediate state, and red the high/soft state. ASM data after MJD 55200 (shown in grey) are affected by 
instrumental decline. The low/hard and intermediate states cannot be distinguished in BAT; data corresponding to these periods of the low/hard or intermediate 
state are therefore shown in black. The lowest panel shows distribution of the available optical spectra from Ondfejov (lower) and Xinglong (upper). The 
spectra are coloured according to the spectral state. 
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Figure 2. Distribution of the optical spectra of HDE 226868 in orbital phase corresponding to different X-ray states. The low/hard state in blue (44 spectra), 
the intermediate state in green (19 spectra), and the high/soft state in red (44 spectra). The longer lines in the upper panel represent a distribution of 17 spectra 
corresponding to the low/hard state in 2003. Similarly in the lower panel, we highlight 20 spectra corresponding to the high/soft state in 2013. 
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Figure 3. Optical spectra showing the orbital phase evolution of the major spectral features - Ha line and He I A6678. The observation dates and orbital 
phases are given in the left and right sides of each spectrum, respectively. The left panel shows the selected spectra corresponding to the low/hard state in 2003. 
The right panel represents the high/soft state in 2013. 


of the result. For the 2003 data, the distribution of spectra across 
the orbital phase is not optimal, showing two large gaps around 
T> = 0.29 and 0.86. However, we still should be able to get rea¬ 
sonable results showing the distribution of emission in the velocity 
space. 

The Doppler tomograms of Ha A6562.76 are shown in 
Figs 4 and 5. The top panels show the input trailed spectra (on the 
left) and reproduced trail from the fitted emission distribution (on 
the right) used to evaluate the goodness of the fit. The four maps 
in the remaining panels correspond, starting top-left and moving 
in clockwise fashion; the recovered Doppler image of the average 
emission from the system, the total modulated emission, the sine 
amplitude map and the cosine amplitude map of the modulated 
emission. The spectra in Fig. 4 from the period 2003-2013 were 
binned to the orbital phase. 20 bins were used in total, with only 
19 of them were filled. The bin 0.4-0.45 had no appropriate data in 
it, hence, it remains white. The Doppler tomogram in Fig. 5 using 
data from 2003 only, was produced with spectra divided in eight 
phase bins - all of them were filled. 

We see that the reconstructed trail in Figs 4 and 5 reproduces 
the original data well in both cases. Poorer orbital phase resolu¬ 


tion in the case the spectra from 2003 is obviously given by the 
highly limited number of spectra available in this period. The basic 
structure of the trail, however, remains clear and comparable to the 
trail constructed from the spectra from the entire observation pe¬ 
riod 2003-2013. The reconstructed Doppler images in both figures 
display, similarly, the same general structure. We notice a conspic¬ 
uous emission spot near the position of the donor star depicted via 
its critical Roche lobe. The two transferring trajectories are indi¬ 
cated along with the positions of the compact companion and the 
centre of the donor which are marked with the ’x’ crosses. The 
barycentre of the system is shown by a cross 

3.2 The high/soft state 

Following are the results of high/soft state data analysis. First, we 
apply the Doppler mapping on the entire ensemble of 44 spectra 
taken in 2003-2013. Then we analyse only the 20 spectra gathered 
in 2013 over a four-month interval. In this case, the situation is 
somewhat better in comparison to the low/hard state in 2013. The 
data, however, also display a couple of large gaps with the most 
prominent being centred around T> = 0.25 and 0.5. 
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Figure 4. Modulation Doppler tomogram of Ha A6562.76 line produced 
using all spectra from the period 2003-2013 appropriate to the low/hard 
state of Cygnus X-1. The spectra were divided into 20 phase bins. The hgure 
shows the trailed spectra on top-left panel, the simulated trail on the top- 
right panel, the distribution of emission in the system in velocity space on 
the middle-left panel, and the variations of the emission on the middle-right 
panel. The bottom panels represent the cosine (on the left) and sine (on the 
right) amplitude maps. 


The results from modulation mapping of Ha line during the 
high/soft state are shown in Figs 6 and 7 for the data covering 2003- 
2013 period and the 2013 data only, respectively. The spectra in 
Fig. 6 were binned to 17 orbital phase bins, with 16 of them filled 
and one remained empty. The smaller number of bins is caused 
with the less ideal orbital phase distribution of the observed data. 
The Doppler tomogram in Fig. 7 uses the 2013 data only, and was 
produced with spectra decided in nine phase bins - all of them were 
filled. 

The reconstructed trails for the high/soft state (Fig. 6 and 7) 
are not as good as for the low/hard state. This is mainly caused by 
the larger gaps in the phase distribution of the data and, by exten¬ 
sion, by the need of having larger phase bins. Another feature wors¬ 
ening the situation is the complex structure of the Ha line trail. This 
makes the analysis significantly more difficult in the high/soft state. 
In general, the Ha emission is less pronounced than in the low/hard 
state. But we still notice a number of bright spots in the emission 
distribution map. These can be remnants of a disc-like structure 
around the donor star which would correspond to the P-Cygni pro¬ 
file coming from the expanding stellar wind. But at least some of 
them are more likely to be numerical artefacts. Despite the lesser 



(km/s) (km/s) 

Figure 5. Modulation Doppler tomogram of Ha A6562.76 line produced 
using spectra from 2003 only appropriate to the low/hard state of Cygnus 
X-1. The spectra were divided into 8 phase bins. Panels as defined in Fig. 4. 

phase coverage and generally poorer quality of the tomograms in 
the high/soft state, we can conclude that the emission representing 
the gas stream from the donor star is no longer there, suggesting 
substantial redistribution of the circumstellar material during the 
transition between individual spectral states. 


4 RADIATION HYDRODYNAMIC MODEL 

While not as intuitive to interpret as spatial distribution in Cartesian 
coordinates, Doppler tomograms in velocity space can be obtained 
without prior assumptions about the velocity field in the fiow as a 
function of position. This significantly simplifies the inversion pro¬ 
cess and allows application of Doppler tomography in a variety of 
conditions where the nature of the fiow is not a-priori known. A dis¬ 
advantage is that two spatially very distinct parcels of gas can have 
precisely the same velocity coordinates; therefore, the transition to 
spatial coordinates is not uniquely defined. 

Aside from this non-uniqueness, it is in principle possible for 
a given model to convert the velocity coordinates into position co¬ 
ordinates. The only thing that is required is a specification of the 
velocity at every point in the system. However, there are two com¬ 
plications. First, the transformation between velocity and position 
is often not known. In fact, it is rarely known, given that it is likely 
that deviations from Keplerian fiow occur. This means that position 
maps would require re-computation each time system parameters 
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Figure 6. Modulation Doppler tomogram of Ha A6562.76 line produced 
using all spectra from the period 2003-2013 appropriate to the high/soft 
state of Cygnus X-1. The spectra were divided into 17 phase bins. Panels as 
defined in Fig. 4. 


Table 2. Simulation parameters 


Parameter 

Symbol 

Used value 

Effective temperature of primary . 

. Teff (K) 

30000 

Average radius of primary. 

. R. (Rq) 

18 

Luminosity of primary . 

L* (Lq) 

2.25 X 10^ 

X-ray luminosity 

Low/hard state . 

^LH(ergs“^) 

1.9 X 10^'^ 

High/soft state. 

^Hs(ergs“^) 

3.3 X 10^'^ 

Mass of donor . 

M* (Mq) 

24 

Mass of BH . 

. Mx (Mq) 

8.7 

Mass-loss rate . 

. M(Moyr-i) 

2 X 10-® 

Orbital period. 

Povh (d) 

5.599829 

Binary separation . 

. D (Rq) 

42.4 

Eccentricity . 

e 

0 


were updated. Secondly, the same place in the system can produce 
emission at more than one velocity. There are many examples of 
bright-spot emission from the gas stream while the disc at the same 
location produces emission at completely different velocity coor¬ 
dinates (cf. e.g. Longa-Pena et al. 2015). If such data are imaged 
into position coordinates on the basis of Keplerian rotation, a spot 


Figure 7. Modulation Doppler tomogram of Ha A6562.76 line produced 
using spectra from 2013 when Cygnus X-1 was in the high/soft state. The 
spectra were divided into nine phase bins. Panels as defined in Fig. 4. 

of emission could be produced at a spurious location in the disc. 
These interpretation difficulties can be potentially resolved only in 
eclipsing binaries. 

While it is not trivial to convert the tomograms into posi¬ 
tion coordinates, there is no difficulty in translating predictions 
of any theoretical model into velocity space. Ideally, the model- 
observation comparison should be made by predicting trailed spec¬ 
tra, abandoning the Doppler maps altogether. However, Doppler 
maps still play a role in those theoretical models that are unable 
to include all peculiarities of real systems, and thus the comparison 
is easier in velocity coordinates. In CH, we presented an enhanced 
radiation hydrodynamic model of the stellar wind in HMXBs first 
introduced by Hadrava & Cechura (2012). We use the model to sim¬ 
ulate the structure and dynamics of the circumstellar environment 
in the vicinity of Cygnus X-1 in both low/hard and high/soft states. 
Then we convert this information into synthetic Doppler maps di¬ 
rectly comparable to the tomograms obtained from the observa¬ 
tions. 

Despite being extensively studied for five decades, key pa¬ 
rameters of Cygnus X-1 remain uncertain. Many estimates have 
been made of the masses of both components of the binary - e.g., 
Caballero-Nieves et al. (2009) inferred wide observational limits 
stating the mass of 23lg Mq and III 3 Mq for the primary and 
the compact component of the system, respectively. However, most 
of these estimates are unreliable because they are based on unsatis¬ 
factory determination of the distance to Cygnus X-1. In their recent 
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Figure 8. Distributions of physical quantities in the equatorial plane of Cygnus X-1 as calculated with the radiation hydrodynamic model of the stellar wind. 
The displayed quantities are averaged over 1 Porb ^^d represent the density distribution, local temperature, and the ratio of neutral to all hydrogen atoms. The 
upper and bottom panels correspond to the low/hard and the high/soft state, respectively. 


paper, Orosz et al. (2011) were able to constrain strongly some of 
the principal parameters of Cygnus X-1. By utilizing an unprece¬ 
dentedly precise distance from a trigonometric parallax measure¬ 
ment of Cygnus X-1 (1.86lo' ii kpc; Reid et al. 2011), they found 
masses of M. = 19.16 ± 1.90 M© and = 14.81 ± 0.98 M© 
for the 0-star and black hole, respectively. These estimates are con¬ 
siderably more robust than previous ones, owing largely to the new 
parallax distance. 

In a recent paper, however, Ziolkowski (2014) showed that the 
mass of the supergiant is inconsistent with the evolutionary mod¬ 
els for massive core hydrogen burning stars. Discrepancies of this 
kind between stellar classification and observed luminosities, in¬ 
ferred masses and temperatures, have previously been reported for 
the donors in HMXBs (cf., e.g., Conti 1978; Kaper 2001; Coe et al. 
2015). Based on the evolutionary models, the mass of the super¬ 
giant is, most likely, in the range of 25 to 35 M©. The correspond¬ 
ing mass of the black hole is in the range of 13 to 23 M©. If, as a 
result of the rotation induced mixing, the hydrogen content of HD 
226868 is equal to about 0.6 (as suggested by some observations), 
then its present mass may be somewhat lower, ~ 24 M©. 

The parameters of the following simulation were chosen in ac¬ 
cordance with the observed parameters of Cygnus X-l/HD 226868 
(see Table 2). Regarding the above summarized uncertainties, we 
chose the component masses as a first probe of their values to con¬ 
vey approximately our data at the present time. We fixed the mass 
of the donor at M* = 24 M© (a value close to the upper limit of the 


conceivable mass range) in order to cover a wider range of the mass 
ratios. Using our optical data, we determined a position of the donor 
in the velocity space by creating a Doppler map of an inverse He I 
absorption line. Because this line originates from the donor itself, 
such an emission will map to a region that represents the velocity 
amplitude of the donor. By adjusting the mass ratio of the system, 
we fitted the size and the shape of the critical Roche lobe on the 
central non-emitting feature in Figs 4 and 5, and thus derived the 
mass of the compact companion equal to = 8.7 M©. 

The orbital period is 5.6 d and the corresponding orbital sep¬ 
aration of the components is 42i^©. We assume the primary to be 
tidally locked and synchronously rotating with the companion or¬ 
bital motion. See CH for the detailed description of the physical 
model we use. 

In Fig. 8, we present the results of our simulations for two 
values of the X-ray luminosity. The upper panels represent the case 
with the X-ray luminosity equal to 1.9 x 10^^ erg correspond¬ 
ing to the low/hard X-ray state of Cygnus X-1. The lower panels 
show the case of the X-ray luminosity equal to 3.3 X lO^^ergs”^ 
which corresponds to the high/soft X-ray state. Each column in 
Fig. 8 represents a different physical quantity. From left to right, 
it is the density distribution p, the local temperature T, and the 
amount of neutral hydrogen nu within the orbital plane of Cygnus 
X-1. The solution of the model in both cases is quasi-stationary. 
While the general structure of the fiow remains steady and is mainly 
determined by the X-ray state, the accreting region which is directly 
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affected by the interaction of the shock with the undisturbed wind is 
rather turbulent and endures quasi-periodic osculations. Since these 
fluctuations are of a random nature and take place on the time- 
scales considerably shorter than the orbital period, we Alter them 
out by time-averaging the values of the quantities of interest over 1 
-Porb* 

The presence of the strong X-ray source has a profound ef¬ 
fect on the structure and dynamics of the wind. Ionizing the gas 
material severely limits the efficiency of the line-driven mecha¬ 
nism introduced by Castor, Abbott and Klein (Castor et al. 1975, 
hereafter CAK), which is the dominant force propelling the wind. 
In the low/hard state, the outflow from the hemisphere facing the 
X-ray source is noticeably slowed down. This is in direct contrast 
with the wind launching from the opposite hemisphere lying in the 
X-ray shadow cast by the donor where the wind is relatively unaf¬ 
fected by the X-ray radiation. The shadow provides an environment 
where the wind can be accelerated without interruption, leading to 
velocities that are an order of magnitude higher than in the slow- 
wind region where the acceleration mechanism is impaired right at 
the base of the wind. Slowing down the wind gives a rise to a large 
accretion disc around the compact object which is preceded by an 
extensive shock. The stellar wind is highly anisotropic. The gas 
launched from the hemisphere facing the X-ray source condensates 
into a narrow dense stream and passes in the vicinity of Li point 
in the direction of the compact companion. The mass-loss rate in 
the direction to the Li point is thus significantly enhanced. In the 
high/soft state, the X-ray radiation forms a bubble of full ionization 
which almost reaches the surface of the donor rendering the CAK 
mechanism ineffective and preventing the wind to achieve the es¬ 
cape velocity. This effectively cuts off the outflow of the material 
from the facing hemisphere. The most prominent consequence of 
the increased X-ray luminosity is the disruption of the narrow dense 
stream between both components of the binary. The mass transfer 
from the donor to the accretion disc is then interrupted and the ac¬ 
cretion disc shrinks as the matter within it is accreted on to the com¬ 
pact companion. Lacking the incoming flow of gas to compress the 
shock enveloping the accretion disc, the shock expands in the wide 
region around the compact companion abruptly increasing the tem¬ 
perature of the remaining gas. Without any prevailing velocity field, 
the movement of gas in this region is mostly slow and turbulent, and 
the bow shock around the accretion disc completely disappears. 


4.1 Production of synthetic Doppler tomograms 

We use the predictions of our model to construct synthetic Doppler 
tomograms mapping the Ha line emission in the low/hard and 
high/soft state of Cygnus X-1. From the density distribution in the 
top left panel of Fig. 8, we expect that the Doppler maps of Cygnus 
X-1 in the low/hard state will show the presence of features cor¬ 
responding to the gas stream and also to the extensive accretion 
disc. Such a disc would manifest itself in the Dopppler map as a 
ring-like structure around the position of the compact companion. 
But we see no signs of the accretion disc present in the spectra of 
HDE 226868 in Figs 4 and 5. In the high/soft state, we expect that 
the increased X-ray emission from the compact companion cuts off 
the gas stream supplying the accretion disc which becomes much 
less pronounced and possibly not visible at all in Ha line. 

To appreciate the conditions within the circumstellar medium 
in different spectral states, it is convenient to display local temper¬ 
ature and the amount of neutral hydrogen nu which we define as a 


ratio 


^ n(H°) 

np+n(H0) 


( 1 ) 


where n(H°) is number density of neutral hydrogen, and np is pro¬ 
ton density. Since an Ha line is emitted by neutral hydrogen, nu 
will serve as a tracer of Ha line emission. 

We can transform the last column in Fig. 8 into velocity space. 
The resultant Doppler maps corresponding to the distribution of 
neutral hydrogen in both major X-ray states are displayed in Fig. 9. 
The colour-scale shows an expected relative flux in Ha line cor¬ 
responding to the distribution of neutral hydrogen in the orbital 
plane. We can now directly compare both panels of Fig. 9 with 
the Doppler tomograms in Figs 4 and 6 produced from the optical 
observations. We see that the occurrence of the highest concen¬ 
tration of neutral hydrogen in the velocity space in the synthetic 
tomogram roughly coincides with the emission region of Ha line. 
As expected, the Ha line emission practically disappears when we 
switched to the high/soft state (right-hand-side panel of Fig. 9). 

In Fig. 8, the hydrodynamic shock around the accretion re¬ 
gion abruptly increases the temperature of the medium up to 10^ 
K. This temperature increase is high enough to cause hydrogen to 
be completely ionized due to mutual collisions. The region around 
the compact companion shows only trace amounts of neutral hy¬ 
drogen leaving the gas stream as the most prominent source of Ha 
line emission in the system. Due to the increased emission of X- 
ray in the high/soft state (lower panels of Fig. 8), the outflow of 
material from the hemisphere facing the compact companion is in¬ 
terrupted. As a consequence, density in the half-space around the 
compact companion is considerably diminished, allowing the hy¬ 
drodynamic shock to expand and heat up the remaining gas. In both 
cases, however (as shown in the rightmost column of Fig. 8), de¬ 
spite the considerable amount of gas there are only traces of neu¬ 
tral hydrogen predicted in the vicinity of the compact companion. 
Thus, the only strong emission of Ha originates in the dense stream 
which is present in the low/hard state. This emission translates to a 
Doppler map as a localized maximum of intensity near the ballistic 
trajectory (left panel of Fig. 9). The fact that the maximum does 
not lie precisely at the ballistic trajectory itself is because such a 
stream would be launched from the Li point at the critical Roche 
equipotential. However, we assume the donor to fill only 95% of 
the critical Roche lobe. Thus the stream is shifted towards the cen¬ 
tre of primary in Vy. The shift towards higher negative Vx is caused 
by the radiative (CAK-) acceleration which is higher compared to 
a free-fall of the ballistic trajectory in the effective Roche poten¬ 
tial. In the case of the high/soft state which is displayed in the right 
panel of Fig. 9, we notice that the maximum of intensity has disap¬ 
peared and the remaining emission is scattered and weak. In either 
case, however, we do not register any contribution of the accretion 
disc to the Ha line emission. 

Besides the expected dominant emission coming from the 
dense matter flow in the low/hard state, there is also a fine struc¬ 
ture of less pronounced ambient emission present in the synthetic 
Doppler map in the left panel of Fig. 9. To get a better idea of 
which parts of the wind contribute to this emission, we attempt to 
associate the emission in the velocity space with its spatial posi¬ 
tion. In Fig. 10, we mark some of the regions of the wind and trace 
them to their corresponding locations in the Doppler tomogram. 
This approach gives us a set of paired points in velocity and spa¬ 
tial coordinates labelled with numbers. The locations of the points 
1 and 2 in velocity space are rather intuitive and coincide with the 
region of enhanced emission. However other points, although quite 
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Low/Hard state 



0 0.25 0.50 0.75 1 



High/Soft state 



0 0.25 0.50 0.75 1 



Figure 9. Synthetic Doppler tomograms of the predicted Ha line emission based on the model data of neutral hydrogen in two major states of Cygnus X-1 (as 
shown in the last column of Fig. 8). The maps are directly comparable with the Doppler maps acquired from the observational data of HDE 226868 in Fig. 4 
and 6. 



Figure 10. Associated regions of the wind in spatial and velocity coordinates labelled with corresponding numbers. Displayed quantities are the distribution 
of neutral hydrogen nn in the orbital plane (the left panel) and corresponding relative Ha line emission in the velocity space (the right panel). 


distant in spacial terms, map into very similar areas (or even overlap 
each other) in the velocity space - e.g. 2-7 or 5-8. This highlights 
the ambiguity of identifying structures in Doppler tomograms with 
their corresponding spatial locations. 


5 CONCLUSION 

In addition to the radiation hydrodynamic simulations of the stel¬ 
lar wind in HMXBs, we also performed optical observations and 
Doppler tomographic analysis of the phase-resolved spectroscopic 
observations of Ha line of Cygnus X-1 in the period of 2003-2013. 
In order to obtain the best orbital phase coverage possible, we put 
together the spectra we observed (in part) using the Ondfejov 2-m 
telescope in the Czech Republic with the spectra taken at the Xin- 


glong station in China. We utilized the X-ray observations from 
RXTE, MAXI and Swift-BAI in order to compartmentalize the 
optical data according to their X-ray spectral state. The results of 
our analysis show dramatic redistribution of the circumstellar ma¬ 
terial in Cygnus X-1 during the transition between the low/hard 
and high/soft states. In the low/hard state, we find an Ha emission 
that is consistent with a dense matter stream (focused stellar wind) 
in the proximity of the inner Lagrangian point Li launching from 
the donor star and accreting on to the compact companion. Despite 
the poorer phase coverage leading to worsened quality of the re¬ 
sulted Doppler maps, we conclude that in the high/soft state when 
the X-ray emission from the compact companion is increased, the 
Ha emission from the matter stream is suppressed. This indicates 
disruption of the mass transfer between the components of the bi¬ 
nary. In the case of the third most prominent X-ray spectral state 
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in Cygnus X-1 - the intermediate state - we did not acquire a suf¬ 
ficient amount of data which would allow us to carry out the to¬ 
mographic analysis. Besides the insufficient phase coverage of the 
optical spectra, we can expect that during this state the distribution 
of the circumstellar matter is erratic and highly variable. In order to 
produce a reliable Doppler tomogram in this X-ray state, we would 
need observations covering no more than a few orbital periods. 

The results of our simulations suggest (see Fig. 4 in CH) that 
the large structures, such as the accretion disc or the hydrodynamic 
shock that envelopes it, are rather insensitive to small variations of 
the mass ratio of the compact component to the donor. As a con¬ 
sequence, the hydrodynamic simulations are not a suitable tool to 
better constrain the system parameters such as masses of the indi¬ 
vidual components. The accretion process, on the other hand, has a 
decisive influence on the circumstellar matter and the structure of 
the wind. The X-ray feedback impacts the ionization structure of 
the wind and thus its ability to be accelerated by the radiation field 
of the donor. There is also evidence for a jet in Cygnus X-1 (Gallo 
et al. 2005) which is currently not included in our simulations but 
which is likely to infiuence the wind dynamics and may introduce 
further shocks. Therefore, we plan to further develop the code to 
enhance its capabilities by adopting an improved physical model. 
For instance, extending the computation of ionization structure of 
the wind material to elements other than hydrogen would be very 
beneficial for our understanding of the wind structure and proper¬ 
ties. 

The jet of the type observed by Gallo et al. (2005) would be 
generated in close proximity to the compact companion. Because 
of the high computational demands, however, we cannot use the 
fine resolution necessary to accurately simulate the physical pro¬ 
cesses in this region. For further advancement in this matter, we 
have to first improve our model of the accretion disc region tak¬ 
ing into account effects of magnetic field and general relativity. 
This general relativistic radiation magneto-hydrodynamic approach 
(GR-MHD) would provide us with a more self-consistent model of 
X-ray binaries. To advance our observation-interpretation method 
which utilizes the calculated synthetic tomograms, we also intend 
to acquire data of sources other than Cygnus X-1. For this pur¬ 
pose, low-mass X-ray binaries with typical hours-long orbital pe¬ 
riods represent ideal candidates because they allow for full orbital 
coverage in a matter of one or a few consecutive nights. 

Our time-dependent simulations show that even though the 
wind medium can be susceptible to turbulent motion at the small- 
scales, the overall solution remains quasi-stationary. This is in 
agreement with our initial assumption that the large-scale structures 
present in the wind are stable, and are mainly determined by the X- 
ray state of Cygnus X-1. Therefore, our use of the optical data that 
span a period much longer than the orbital period of Cygnus X- 
1 for Doppler tomography, seems to be well founded. We use the 
results of our model of Cygnus X-1 to calculate the amount of neu¬ 
tral hydrogen in the circumstellar matter which is responsible for 
generating the Ha line emission. For this distribution, we produce 
synthetic Doppler tomograms corresponding to both major X-ray 
states - the low/hard and high/soft state - and compare them with 
the tomograms obtained from the spectroscopic analysis. The com¬ 
parison of the observational data with the numerical simulations is 
done in velocity rather than spatial coordinates to avoid the prob¬ 
lems with ambiguity. In both X-ray states, we find the predicted 
distributions of neutral hydrogen to be in good agreement with the 
observations, thus confirming the assumed parameters of Cygnus 
X-1. 
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APPENDIX A: 

The complete list of the spectra used in the Doppler analysis is summarized in Table Al. 


Table Al: The list of spectra used for the Doppler tomography 


Date 

(yyyy/mm/dd) 

UT 

start 

(hh:mm:ss) 

Exposure 

time 

(s) 

Heliocentric 

Julian date 

Wavelength 

range 

(A) 

Orbital 

phase 

Spectral 

state* 

2003/04/01 

03:27:09 

2000 

2452730.6420 

6257-6769 

0.613 

LH 

2003/04/17 

02:14:28 

3466 

2452746.5924 

6258-6770 

0.463 

LH 

2003/07/14 

21:55:29 

4000 

2452835.4167 

6258-6771 

0.326 

HS 

2003/07/20 

23:39:47 

3600 

2452841.4893 

6258-6771 

0.231 

HS 

2003/07/25 

23:16:59 

3600 

2452846.4735 

6258-6771 

0.300 

HS 

2003/08/05 

22:28:34 

600 

2452857.4399 

6258-6771 

0.255 

I 

2003/08/05 

22:43:14 

3600 

2452857.4501 

6258-6771 

0.260 

I 

2003/08/08 

23:54:02 

1800 

2452860.4993 

6259-6771 

0.624 

I 

2003/08/09 

00:30:30 

1800 

2452860.5246 

6259-6771 

0.807 

I 

2003/08/09 

01:04:10 

1800 

2452860.5480 

6259-6771 

0.811 

I 

2003/08/10 

00:21:41 

3600 

2452861.5185 

6258-6771 

0.987 

LH 

2003/08/11 

01:01:03 

300 

2452862.5458 

6259-6771 

0.167 

LH 

2003/08/11 

01:08:55 

3600 

2452862.5513 

6259-6771 

0.171 

LH 

2003/08/22 

22:59:45 

1800 

2452874.4614 

6259-6772 

0.296 

I 

2003/08/22 

23:32:08 

1800 

2452874.4839 

6259-6772 

0.300 

I 

2003/08/23 

00:03:43 

1800 

2452874.5059 

6259-6772 

0.304 

I 

2003/08/23 

00:34:51 

1800 

2452874.5275 

6259-6772 

0.308 

I 

2003/08/24 

20:20:45 

1800 

2452876.3510 

6259-6772 

0.632 

I 

2003/08/24 

20:54:48 

1800 

2452876.3746 

6259-6772 

0.636 

I 

2003/08/24 

21:27:09 

1800 

2452876.3971 

6259-6772 

0.640 

I 

2003/08/24 

22:01:29 

1800 

2452876.4209 

6259-6772 

0.644 

I 

2003/08/25 

01:26:12 

1800 

2452876.5631 

6259-6772 

0.669 

I 

2003/08/25 

02:04:27 

1800 

2452876.5897 

6259-6772 

0.674 

I 

2003/08/26 

21:52:08 

400 

2452878.4144 

6259-6772 

0.000 

LH 

2003/08/26 

22:02:47 

3300 

2452878.4218 

6259-6772 

0.005 

LH 

2003/09/15 

22:33:24 

4000 

2452898.4425 

6258-6770 

0.581 

LH 

2003/09/16 

22:20:30 

3600 

2452899.4335 

6258-6770 

0.754 

LH 

2003/09/18 

19:57:43 

4000 

2452901.3343 

6261-6773 

0.097 

LH 

2003/09/20 

22:13:42 

1200 

2452903.4286 

6259-6772 

0.468 

LH 

2003/09/20 

22:37:49 

5175 

2452903.4454 

6259-6772 

0.475 

LH 

2003/09/21 

22:41:52 

5499 

2452904.4481 

6259-6771 

0.655 

LH 

2003/09/24 

23:19:04 

300 

2452907.4756 

6258-6770 

0.190 

LH 

2003/09/24 

23:26:23 

3600 

2452907.4998 

6258-6770 

0.194 

LH 

2003/09/25 

23:44:19 

180 

2452908.4924 

6258-6770 

0.371 

LH 

2003/09/25 

23:49:54 

3077 

2452908.5130 

6258-6770 

0.375 

LH 

2004/04/12 

00:33:19 

7200 

2453107.5636 

6257-6769 

0.921 

LH 

2004/08/09 

00:20:20 

6300 

2453226.5540 

6264-6776 

0.170 

LH 

2004/08/24 

22:46:10 

300 

2453242.4537 

6263-6776 

0.009 

LH 

2004/08/24 

22:54:29 

3600 

2453242.4786 

6263-6776 

0.013 

LH 

2004/09/02 

20:13:36 

300 

2453251.3475 

6263-6776 

0.597 

LH 

2004/09/02 

20:24:29 

1800 

2453251.3638 

6263-6776 

0.600 

LH 

2004/09/21** 

14:46:04 

900 

2453270.1205 

6338-6718 

0.949 

HS 

2004/09/21** 

15:01:51 

900 

2453270.1315 

6338-6718 

0.951 

HS 

2004/09/22** 

14:36:11 

900 

2453271.1137 

6337-6811 

0.126 

HS 

2004/09/25** 

13:47:11 

1000 

2453274.0802 

6338-6738 

0.656 

HS 

2004/09/25** 

14:04:44 

1000 

2453274.0924 

6338-6738 

0.658 

HS 

2004/09/26** 

13:23:08 

1000 

2453275.0635 

6340-6700 

0.831 

HS 

2004/09/26** 

13:40:42 

1000 

2453275.0757 

6340-6700 

0.834 

HS 

2004/10/24 

19:01:57 

300 

2453303.2954 

6266-6779 

0.874 

I 

2004/10/24 

19:12:04 

1800 

2453303.3112 

6266-6779 

0.877 

I 

2004/11/24 

16:50:00 

300 

2453334.2020 

6264-6776 

0.393 

I 

2004/11/24 

16:57:04 

2400 

2453334.2191 

6264-6776 

0.396 

I 

2006/09/26** 

13:55:43 

1200 

2454005.0873 

6339-6777 

0.196 

LH 

2006/09/26** 

14:17:18 

1200 

2454005.1023 

6339-6777 

0.199 

LH 

2006/09/27** 

13:16:20 

1200 

2454006.0600 

6339-6777 

0.370 

LH 

2006/09/27** 

13:55:07 

1200 

2454006.0869 

6339-6777 

0.375 

LH 

2006/09/28** 

12:14:02 

1200 

2454007.0167 

6339-6866 

0.541 

LH 

2006/09/28** 

13:15:46 

1200 

2454007.0596 

6340-6867 

0.549 

LH 
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Date 

(yyyy/mm/dd) 

UT 

start 

(hh:mm:ss) 

Exposure 

time 

(s) 

Heliocentric 

Julian date 

Wavelength 

range 

(A) 

Orbital 

phase 

Spectral 

state* 

2006/09/29** 

12:19:54 

1200 

2454008.0138 

6338-6787 

0.720 

LH 

2006/10/01** 

11:59:59 

600 

2454010.0035 

6335-6798 

0.075 

LH 

2006/10/01** 

15:04:09 

600 

2454010.1314 

6334-6797 

0.098 

LH 

2006/10/02** 

12:04:59 

600 

2454011.0069 

6333-6797 

0.254 

LH 

2006/10/02** 

14:30:00 

600 

2454011.1076 

6332-6796 

0.272 

LH 

2007/07/06 

22:49:27 

3600 

2454288.4748 

6257-6769 

0.804 

LH 

2008/04/24 

21:24:46 

7200 

2454581.4334 

6251-6764 

0.120 

LH 

2008/05/06 

20:53:11 

7200 

2454593.4121 

6251-6764 

0.259 

LH 

2008/05/07 

00:28:20 

7200 

2454593.5615 

6251-6764 

0.286 

LH 

2008/06/02 

23:48:45 

7200 

2454620.5355 

6249-6762 

0.103 

LH 

2008/09/09 

21:36:53 

8000 

2454719.4497 

6252-6765 

0.766 

LH 

2009/08/14 

22:07:44 

3600 

2455058.4462 

6253-6765 

0.303 

LH 

2010/04/27 

00:53:04 

2600 

2455313.5515 

6261-6774 

0.859 

LH 

2010/09/17 

20:53:54 

4200 

2455457.3976 

6256-6768 

0.547 

HS 

2010/10/08 

19:07:35 

3600 

2455478.3194 

6255-6767 

0.283 

HS 

2011/05/05 

20:52:34 

1836 

2455687.3805 

6252-6765 

0.616 

LH 

2011/10/21 

18:45:58 

1800 

2455856.2932 

6252-6764 

0.780 

HS 

2011/10/21 

19:19:18 

1800 

2455856.3164 

6252-6764 

0.784 

HS 

2011/10/21 

21:30:30 

1800 

2455856.4075 

6252-6764 

0.801 

HS 

2011/10/21 

22:03:19 

1800 

2455856.4303 

6252-6764 

0.805 

HS 

2011/10/22 

21:59:24 

1800 

2455857.4275 

6251-6763 

0.983 

HS 

2011/10/22 

22:32:19 

1500 

2455857.4486 

6251-6763 

0.987 

HS 

2011/11/12 

16:25:29 

1800 

2455878.1944 

6252-6764 

0.691 

HS 

2011/11/12 

19:06:26 

1800 

2455878.3062 

6252-6764 

0.711 

HS 

2011/11/13 

17:51:39 

1800 

2455879.2542 

6252-6765 

0.881 

HS 

2011/11/13 

19:40:10 

1800 

2455879.3295 

6252-6765 

0.894 

HS 

2011/12/10 

16:56:27 

1800 

2455906.2144 

6252-6765 

0.695 

LH 

2012/08/20 

22:21:20 

1031 

2456160.4408 

6253-6766 

0.094 

HS 

2012/08/27 

18:58:18 

4500 

2456167.3197 

6253-6766 

0.322 

HS 

2013/07/17 

00:47:50 

1801 

2456490.5469 

6263-6736 

0.043 

HS 

2013/07/18 

20:52:57 

2001 

2456492.3850 

6263-6736 

0.372 

HS 

2013/07/19 

01:55:16 

1394 

2456492.5914 

6263-6736 

0.408 

HS 

2013/08/16 

22:35:10 

3001 

2456521.4618 

6262-6735 

0.564 

HS 

2013/08/17 

23:15:07 

3001 

2456522.4895 

6262-6735 

0.748 

HS 

2013/08/22 

23:48:31 

3001 

2456527.5127 

6262-6735 

0.645 

HS 

2013/08/23 

21:20:18 

6203 

2456528.4283 

6262-6735 

0.808 

HS 

2013/08/24 

19:00:15 

5256 

2456529.3255 

6262-6734 

0.968 

HS 

2013/08/28 

22:16:02 

2307 

2456533.4443 

6262-6734 

0.704 

HS 

2013/08/29 

22:53:50 

3001 

2456534.4746 

6262-6734 

0.888 

HS 

2013/09/04 

23:07:07 

3001 

2456540.4836 

6262-6734 

0.961 

HS 

2013/09/09 

20:52:24 

2701 

2456545.3882 

6262-6734 

0.837 

HS 

2013/09/27 

20:17:04 

8991 

2456563.3993 

6262-6735 

0.053 

HS 

2013/09/29 

21:34:03 

2501 

2456565.4151 

6262-6735 

0.413 

HS 

2013/09/30 

21:40:51 

3601 

2456566.4262 

6262-6734 

0.594 

HS 

2013/10/06 

20:30:59 

6282 

2456572.3929 

6262-6734 

0.659 

HS 

2013/10/07 

22:23:27 

5114 

2456573.4642 

6262-6734 

0.850 

HS 

2013/10/08 

18:55:12 

2968 

2456574.3071 

6262-6735 

0.001 

HS 

2013/10/19 

20:23:17 

4037 

2456585.3739 

6262-6735 

0.977 

HS 

2013/10/31 

20:43:25 

3795 

2456597.3858 

6262-6735 

0.122 

HS 


* low/hard state (LH), high/soft state (HS) and intermediate state (I) 
** spectra observed at the Xinglong station (China) 



